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ABSTRACT 


An acoustic prediction capability for supersonic axisymmetric jets was developed on the basis of 
OVERFLOW Navier-Stokes CFD (Computational Fluid Dynamics) code of NASA Langley Re- 
search Center. Reynolds-averaged turbulent stresses in the flow field are modeled with the aid of 
Spalart-Allmaras one-equation turbulence model. Appropriate acoustic and outflow boundary 
conditions were implemented to compute time-dependent acoustic pressure in the nonlinear 
source-field. Based on the specification of acoustic pressure, its temporal and normal derivatives 
on the Kirchhoff surface, the near-field and the far-field sound pressure levels are computed via 
Kirchhoff surface integral, with the Kirchhoff surface chosen to enclose the nonlinear sound 
source region described by the CFD code. The methods are validated by a comparison of the 
predictions of sound pressure levels with the available data (Seiner et al., 1992) for an axisym- 
metric turbulent supersonic (Mach 2) perfectly expanded jet. 
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ABBREVIATIONS, ACRONYMS, AND SYMBOLS 


CAA 

Computational Aeroacoustics 

CFD 

Computational Fluid Dynamics 

dB 

decibel 

FFT 

Fast Fourier Transform 

gn 2 

Gaseous nitrogen 

Hz 

hertz 

kHz 

kilohertz 

K 

kelvin 

KSC 

Kennedy Space Center 

m 

meter 

NASA 

National Aeronautics and Space Administration 

OASPL 

overall sound pressure level (ref. 20 p Pa) 

s 

second 

SPL 

sound pressure level (ref. 20 p Pa) 


Nomenclature 

c=sound speed, -yjyg c R u T IW , ft/s 

dj = nozzle exit diameter, ft 

F = total thrust, lbf 
/ = frequency, Hz 

g r = gravitational constant, (lbm-ft/lbf.s 2 ) 

I = specific impulse (sec) 

M = Mach number, Vic 
m = instability mode 
p = static pressure, lbf/ft 2 

2 

p'= acoustic pressure disturbance (instantaneous pressure), p- p m , lbf/ft 

Q = vector of disturbance variables 

Re = Reynolds number, pu j d J / p 

R u = universal gas constant, lbf.ft/(lbm-mole R) 

r- radial distance 

Tj = nozzle exit radius, ft 

St = Strouhal number, fd J tu j 

T = temperature, R 
V = velocity, m/s 
u = axial velocity, ft/s 
u c = jet center-line velocity at exit, ft/s 
w= radial velocity, ft/s 


IV 



W = molecular weight, lbm/lbm-mole 
x= axial distance from the nozzle exit, ft 
z = radial distance from jet axis, ft 


Greek Symbols 

p = dynamic viscosity, Ibm.ft/s 
p - density, lbm/ft 
y = isentropic exponent 
O) = circular frequency, radians/s 
if/ = angle from jet axis, deg 

Subscripts 

c - centerline 
e = nozzle exit 
7= jet 

t = stagnation or total, or nozzle throat 
oo = ambient 

Superscripts 


prime 


disturbance from the mean value 




1. INTRODUCTION 


Acoustic loads in a launch vehicle environment represent a principal source for inducing struc- 
tural vibration and may be critical to the proper functioning of vehicle components and ground 
support structures and equipment in the immediate vicinity of the launch pad. A knowledge of 
acoustic loads, including the overall sound pressure level (OASPL), sound pressure level (SPL) 
spectrum, and the distribution (or correlation) of surface acoustic loads, is necessary to provide 
the input for vibroacoustic analysis and evaluation of structural integrity. In the design of launch 
vehicles, it is highly desirable that data on acoustic loads (near-field and far-field noise levels) be 
generated both analytically and from testing of small-scale and full-scale models. Since full- 
scale acoustic and vibration testing is often cost prohibitive, the option of small scale testing 
combined with analysis methods remains as a practical alternative. Accurate characterization of 
acoustic loads on launch pad structures thus proves to be a formidable challenge. 

Noise from subsonic jets is mainly due to turbulent mixing, comprising the contributions of 
large-scale and fine-scale structures (Lighthill, 1952, 1954). The turbulent mixing noise is 
mainly broadband. In perfectly expanded supersonic jets (nozzle exit plane pressure equals the 
ambient pressure), the large-scale mixing noise manifests itself primarily as Mach wave radia- 
tion, caused by the supersonic convection of turbulent eddies with respect to the ambient fluid. 

In imperfectly expanded supersonic jets, additional noise is generated on account of broadband 
shock noise and screech tones. 

The prediction of noise generation and propagation from exhaust jets by computational methods 
entails certain numerical requirements (Tam 1995). Although the existing CFD codes can pre- 
dict the noise generation in the nonlinear source-field, their application to acoustic predictions 
for obtaining near-field and far-field sound levels results in inaccuracies due to numerical dis- 
persion and dissipation over long propagation distances. Specialized Computational Aeroacous- 
tics (CAA) methods were developed to resolve these numerical issues satisfactorily (Tam and 
Webb 1993, Lele 1992). These computational methods for direct calculation of noise in the 
near- and the far-field, however, require prohibitively large computational resources. 

This situation motivated in recent years the application of CFD codes to model the source region 
(containing the sound sources), with the near- and far-field sound predicted by either Lighthill's 
acoustic analogy (Lighthill 1952, 1954) or the Kirchhoff surface integral (Kirchhoff 1883). In 
Lighthill's acoustic analogy the Navier-Stokes equations are rearranged in wave-like form with 
Lighthill's stress tensor acting as the source term, and the near-field and far-field sound is ob- 
tained via a volume integral. This volume integral again requires large computational resources. 
On the other hand, the Kirchhoff method requires only a surface integral, with the Kirchhoff 
surface enclosing the nonlinear source region, and is chosen in a region where the linear wave 
equation is valid. Additionally, the Kirchhoff method does not suffer dissipation and dispersion 
errors when the near-field and far-field sounds are directly calculated with the CFD codes. By 
restricting the use of CFD methods to the nonlinear near-field region for source identification 
and employing Kirchhoff integration for the linear region, computational requirements are 
greatly reduced. 
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The ability to accurately predict the turbulent mixing noise sources remains as a key require- 
ment in the overall acoustic analysis. Both the large-scale turbulence structures (instability 
waves) and the fine-scale turbulence represent the mixing noise sources. The large turbulence 
structures are responsible for Mach wave emission that is directional, and the fine-scale turbu- 
lence mixing noise is fairly uniform independent of direction. In the past the two-equation 
k-e turbulence model is generally used to model the sound sources. However, this model en- 
tails large computational requirements (two additional equations) and is often plagued with con- 
vergence problems in flows with shock waves and strong gradients. Furthermore, the universal 
constants in the turbulence model are often altered to fit the theory with the acoustic measure- 
ments. On the other hand, the simple algebraic mixing length type turbulence model cannot ac- 
curately describe the growth of the mixing layers. Thus a one-equation turbulence model seems 
to be an attractive choice for engineering applications. This model generally results in faster 
convergence and requires less computational resources. 

The objective of this investigation is to predict supersonic jet acoustics on the basis of CFD/ 
Kirchhoff formulation in conjunction with a one-equation turbulence model. This report summa- 
rizes the development of a computational procedure in which the OVERFLOW CFD Navier- 
Stokes code (Buning et al. 1998) predicts the nonlinear source-field within the jet and application 
of Kirchhoff surface integral for acoustic radiation. 

2. ANALYSIS 

This method of analysis is based on the application of a CFD code for identifying the noise 
sources in the source-field and Kirchhoff surface integral for the propagation of sound radiation 
to the near-field and the far-field. 

2.1 CFD ANALYSIS 

2.1.1 NUMERICAL ALGORITHM. For the CFD analysis, the OVERFLOW code (Buning et 
al., 1998) of NASA Langley is used for computing the instantaneous flowfield. This code solves 
the three-dimensional compressible turbulent flow Navier-Stokes equations in generalized coor- 
dinates. The one-equation turbulence model of Spalart-Allmaras (1992) is considered. 

2.1.2 BOUNDARY CONDITIONS. To simulate the exterior (infinite boundary) acoustic prob- 
lem, we need to prescribe outflow and radiation boundary conditions at the boundaries of a finite 
computational domain. These boundary conditions must be sufficiently transparent to the outgo- 
ing disturbances so that they exit the computational domain without significant (nonphysical) 
reflections. It is assumed that the boundaries are sufficiently far from the sources. 

The OVERFLOW code, in its present form, is designed primarily for the prediction of steady or 
unsteady turbulent flowfields and used widely to model aerodynamic flows. It does not however 
provide the necessary time-dependent boundary conditions for handling the reflection-free 
acoustic propagation. Therefore appropriate modifications to the code were made to provide a 
set of time-dependent reflection-free boundary conditions (includes periodic inflow, outflow, and 
radiation). 
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a. Inflow Boundary Conditions . At the supersonic inflow, all data are specified, 
such that all time variations of the characteristic variables are set to zero because 
all waves are incoming waves. 

b. Outflow Boundary Conditions . In general at the outflow boundary, the outgoing 
disturbances consist of acoustic, entropy, and vorticity waves. Here we follow 
Thompson’s approach (Thompson 1990) based on one-dimensional characteris- 
tics analysis (Riemann variables). The amplitude of the outward propagating 
waves at the boundary are specified entirely from the solution within the compu- 
tational domain, while the amplitude of the inward propagating waves is pre- 
scribed at the boundary. For details on the method, the original reference should 
be consulted. 

c. Radiation Boundary Condition . At boundaries where there exist only outgoing 
acoustic waves, Tam and Webb (1993) developed a set of radiation boundary 
conditions based on an asymptotic analysis (valid for large distances). These 
boundary conditions, originally formulated in a spherical coordinate system, when 
applied to a cylindrical coordinate system, are given by (Hixon et al., 1995) 


where 


q : = -v 


R R R 


Q' = [p',u',w',p'Y , V = 




(la) 


(lb) 


In the above equations, Q' is the vector of disturbance variables, x the axial distance from the 
origin, and rthe radial distance from the jet axis. The subscripts on Q' indicate partial deriva- 
tives. 


2.2 KIRCHHOFF CODE 


For the acoustic radiation, the Kirchhoff code YORICK developed by Pilon and Lyrintzis (1998) 
of Purdue University was considered. The Kirchhoff surface, enclosing the nonlinear source re- 
gion, is chosen in a region where the linear wave equation is valid. Outside the Kirchhoff sur- 
face, the flow is governed by the linear equation 



a a V 

dt dx J 


P 


= 0 
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where u x is the freestream (ambient) velocity. In the special case of stationary ambient 
( M „> = 0 ), ^ reduces to the simple wave equation. A solution to the pressure field can be ex- 
pressed by the surface integrals as (Morino and Tseng 1990, Lyrintzis & Mankbadi 1996) 


p'(x,y,z,t) = -^~ 

An 


P_ d r o | 1 dp | 1 

r o dn 0 r 0 dn 0 c„r 0 /3 2 dt 


OUl-m dx ° 

dn n 


\ 


dn 0 


dS n 


) 


(2b) 


where 


_ _ [r 0 -M,(x-x')] 
‘ c„0 2 ’ 


/) 2 + (z-zTF 

*4 -Mir 


In the above equations, is the freestream Mach number, c_ the speed of sound in the 
freestream, the prime a point on the Kirchhoff surface, t the retarded time T = t-t' , n 0 the out- 
ward normal to the Kirchhoff surface S 0 . Thus the pressure at any instant in the region outside 

the Kirchhoff surface can be expressed in terms of the information prescribed on the Kirchhoff 
surface. The required data on pressure, and its normal and temporal derivatives on the Kirchhoff 
surface are taken from the CFD solution. 


3. VALIDATION TEST CASE 

The CFD/Kirchhoff methodology formulated above is tested by comparing the predictions with 
the well-known acoustic test data of Seiner et al. (1992) for Mach 2 perfectly expanded super- 
sonic jets of air. The nozzle exit radius is 0.04572 m. Although the data cover a wide range of 
jet total temperature (using propane burners), detailed validation is performed here for a jet total 
temperature of 755 kelvin (K). Validation at other jet total temperatures will be reported in a 
later report. Table 1 shows the data at the nozzle exit. The jet exit velocity is 822 meters per 
second (m/s), and exhausts into quiescent surroundings at a temperature of 280 K. Figure 1 
shows the location of the far-field measurements of the OASPL with which comparisons with the 
CFD results are reported. 

4. NUMERICAL SOLUTION 

4.1 CFD SOLUTION 

4.1.1 AXISYMMETRIC GRID. For the CFD computations, an axisymmetric grid of size 200 x 
100 x 3 (200 grid points in the axial direction, 100 grid points in the radial direction, and 3 in the 
azimuthal plane) is considered (figure 2). A grid length of 60 jet radii and a grid radius of 10 
radii are considered. The grid is clustered in the radial direction near the nozzle wall to resolve 
the shear layer. Because of the perfectly expanded condition, uniform grid is considered ade- 
quate in the axial direction. 
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4.1.2 FLOWFEELD SOLUTION. At the inlet, the flow variables are specified based on the 
nozzle exit conditions. In addition, a periodic disturbance (axisymmetric mode) of a single fre- 
quency is also prescribed as follows: 


u(z) = Asin(tur)exp 



(3) 


where a) is the circular frequency of the disturbance and A the amplitude of disturbance. Based 
on the work of Mankbadi et al. (1998), values of b - 0.1 and h = 0.78361 are considered. The 
amplitude A is taken to be of the order of 0.01. A value of St = 0.2 (based on typical peak SPL 
frequency for supersonic jets) is considered here for obtaining the disturbance frequency a> . 

Appropriate time-dependent boundary conditions were applied to ensure reflection-free bounda- 
ries. An outflow boundary condition of Thompson type (Thompson 1990) is applied at the out- 
flow boundary, which maintains the mean static pressure at the ambient value. An acoustic ra- 
diation boundary condition of Tam and Webb (1993) is applied at the lateral boundary. 

The CFD solution converges after about 8,000 time-step iterations (the code is run in a time- 
accurate manner) before a periodic state is established. An initial steady-state solution (based on 
local time-stepping, requiring about 600 iterations for convergence) has served as the starting 
solution for the time-accurate run. The solution is obtained on an IRIX workstation (SGI Indigo 
machine). About 15 hours of computing time are required for the solution to achieve a periodic 
state. 

4.2 ACOUSTIC SOLUTION 

After a periodic state has been established, the appropriate data from CFD, corresponding to a 
cycle, is communicated to the Kirchhoff code. The Kirchhoff surface is a cylindrical surface co- 
inciding with a grid line. The radius of the Kirchhoff surface is taken as about 6 radii from the 
jet axis. Only the lateral surface of the cylinder is taken into account, and the cylinder ends are 
neglected due to the effects of nonlinearities. The data to be specified on the Kirchhoff surface 
include the instantaneous pressure, the pressure-time derivative, and the pressure gradient normal 
to the surface. The Kirchhoff code then computes the OASPL and the pressure-time signals in 
the near- and the far-field. From the pressure-time signals, it is possible to compute the SPL 
spectrum at any location with the aid of the Fast Fourier Transform (FFT). 

5. RESULTS AND COMPARISON 

5.1 JET MEAN FLO WFIELD 

Figure 3 shows the computed mean velocity contours in the jet. The potential core and the mix- 
ing of the turbulent shear layer are clearly displayed. A comparison of the predicted centerline 
mean axial velocity with the data is presented in figure 4. Both the data and the CFD solution 
suggest that the jet potential core extends to about 15 jet radii from the nozzle exit plane. The jet 
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potential core may be defined as the axial distance where the jet centerline velocity remains at 
about 99% of the jet exhaust velocity. It is known that the major noise production region occurs 
near the end of the potential core, where the most highly amplified instability wave reaches its 
maximum growth (Seiner et al., 1992). 

In the velocity decay region past the core, the CFD solution predicts a slower growth relative to 
the measurements. This discrepancy is attributable to the accuracy of the Spalart-Allmaras one- 
equation turbulence model. Also, the effects of compressibility are not appropriately taken into 
account in this turbulence model, as it was originally developed for incompressible flows. Data 
on compressible shear layers indicate that compressibility tends to reduce the growth rate of the 
shear layers (Goebel and Dutton, 1991) relative to the incompressible case, with an increase in 
the jet convective Mach number. For example, before the potential core closes, a jet convective 
Mach number M c , which may be defined by (Freund, Lele and Moin, 2000) 

M = 7 

The reduced transverse turbulence intensity is believed to be responsible for the observed reduc- 
tion in the spreading rate for compressible mixing layer. By taking into account the effects of 
compressibility on the turbulent kinetic energy and its specific dissipation rate into a two- 
equation turbulence model, the reduction in the mixing length growth rate is found to agree with 
the data (Kandula and Wilcox 1995). 

Figures 5 and 6 respectively display the mean static pressure and mean Mach number along the 
centerline of the jet. The pressure oscillation along the centerline is relatively small and shows 
the accuracy of the boundary conditions. The centerline Mach number variation follows the 
trend of the centerline velocity, as is to be expected. 

Calculations show that with an increase in jet temperature, the potential core length decreases. 
These results are in qualitative agreement with the data of Seiner et al. (1992). As the jet tem- 
perature increases, the convective Mach number increases, and the Mach wave radiation pattern 
is modified, with the OASPL values generally increasing with an increase in jet temperature. 

5.2 JET INSTANTANEOUS PRESSURE AND OASPL CONTOURS 

Figure 7 displays the instantaneous pressure contours in the source-field, showing the directivity 
of Mach wave emission. The lack of significant reflection near the boundaries indicates that the 
outflow and the radiation boundary conditions implemented in the CFD code are satisfactory. A 
display of the OASPL within the jet, characterizing the nonlinear source-field, as directly com- 
puted by the CFD code, is presented in figure 8. Inside the jet shear layer, the pressure levels are 
high (as much as 180 dB), with the highest OASPL occurring about 7 jet radii downstream of the 
exit. 
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5.3 NEAR- AND FAR-FIELD ACOUSTIC SOLUTIONS 


Figure 9 presents the OASPL contours in the near- and the far-field as computed by the 
CFD/Kirchhoff formulation. It is seen that the OASPL ranges from 80 dB to 130 dB in this re- 
gion. The directivity of the Mach wave radiation (emission) is evident. It is generally accepted 
that the mechanism for Mach wave radiation is strongly connected with the amplification of in- 
stability modes, and this amplification scales with the jet convective Mach number. 

The variation of the OASPL in the streamwise direction at z/r y = 80is compared in figure 10 

with the measurements of Seiner et al. (1992). The predicted peak angle of emission of about 
128 degrees is in good agreement with the data, which indicates a peak at 127 degrees. Also, the 
peak value of the OASPL is about 7 dB less than the measured value. There is an underpredic- 
tion of the OASPL at the off-peak locations both upstream and downstream of the peak. This 
underprediction is primarily due to the axisymmetric simulations (axisymmetric disturbance 
mode m = 0 instability) of the present investigation. It is impossible to develop a three- 
dimensional helical mode instability with axisymmetric simulations. Measurements show that 
near the nozzle exit, axisymmetric modes seem to be dominant, and farther downstream the heli- 
cal nature of the large-scale structure is dominant (Viswanathan and Sankar 1995). Also, both 
experiments and hydrodynamic stability analysis suggest that the m = +1 helical mode is domi- 
nant above a jet Mach number of about 1.4. Thus by including azimuthal mode (m = +1) insta- 
bility in a three-dimensional analysis, predictions of sound pressure levels could be considerably 
improved (Garnet and Estivalezes 1998). Such three-dimensional simulations require considera- 
bly more computing resources, which are outside the scope of the present work. Further im- 
provements in the predictions can be realized by considering inflow disturbances of various fre- 
quencies (or random disturbances). 

Figure 1 1 shows the corresponding comparison with the OASPL plotted as a function of the ax- 
ial distance from the exit plane. The computed results of OASPL for three values of z/r y =30, 

55 and 80 are displayed in figure 12, the last one corresponding to the data by Seiner (1992). 
With increasing values of zlr j , the peak value of OASPL decreases, and the corresponding peak 

angle decreases. This is to be expected in view of the directivity of Mach wave radiation. The 
correctness of this trend also serves to verify the Kirchhoff code implementation. 

The spectral content of SPL (dB/Hz) along a direction corresponding to the maximum radiation 
( xtr ] = 60, z / rj = 80 ) is displayed in figure 13. Figure 14 indicates the variation of 1/3 octave 

SPL (dB) as a function of the frequency. A well-defined peak SPL value of about 128 dB is seen 
to occur at a frequency of about 2 kilohertz (kHz), which corresponds to a Strouhal number of 
0.2. This value is consistent with the frequency of the imposed disturbance at the nozzle exit. 
Figure 15 presents the result of figure 14 in a semi-log plot. 
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6. CONCLUSIONS 


The CFD/Kirchhoff approach is capable of predicting the SPL spectrum and the OASPL at any 
location in the near-field and the far-field. The relative reduction in the computational time due 
to the Kirchhoff method and the axisymmetric simulation make it an attractive engineering tool. 
The one-equation turbulence model of Spalart-Allmaras for computing turbulent stresses affords 
additional reduction in the computational time and at the same time improves the robustness of 
the code relative to two-equation turbulence models. The peak angle of Mach wave radiation 
and the peak level of OASPL are satisfactorily predicted. Predictions at off-peak conditions 
could be substantially improved by considering three-dimensional disturbances (azimuthal mode 
instability) and disturbances of various frequencies. 
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Table 1. Summary of Nozzle Conditions 


Stagnation temperature, K 

755 

Nozzle exit diameter, m 

0.09144 

Exit pressure, atm 

1.0 

Exit temperature, K 

419 

Exit velocity, m/s 

822 

Acoustic velocity at exit, m/s 

411 

Nozzle exit Mach number 

2.0 


kesemb 

Ambient pressure, atm 

1.0 

Ambient temperature, K 

280 
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Figure 1. Schematic of the Test Case Nozzle 



Figure 2. Axisymmetric grid (200x 100) 
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Figure 8. Contours of OASPL in the Source-Field as Computed Directly From CFD 
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Figure 9. Contours of OASPL Predicted From CFD/Kirchhoff Method 



Figure 10. Comparison of Angular OASPL Distribution at z ! r =80 
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Figure 11. Comparison of Axial OASPL Distribution at z!r ) - 80 



Figure 12. Computed Axial Distributions of OASPL for z / =30, 55, and 80 




SPL (dB), 1/3 octave 



Figure 13. Computed SPL Distribution at x!r j =60, and zlr } = 8C 



Figure 14. Computed SPL Distribution at x/r =60, and z/r =80 
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Figure 15. Computed SPL Distribution (Semi-Log Plot) at x/r, =60, and z/r 
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